Quantum Phases of Ultracold Bosonic Atoms in two Bands of an Optical-Lattice 

coupled by a Cavity Field 
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We study the quantum phase transitions between superfluid and Mott insulator states for ultracold 
bosons occupying two bands of an optical lattice. The two atomic states are resonantly coupled by a 
single cavity mode which mediates transitions between the two bosonic particle modes via absorption 
or emission of a cavity photon. This coupling between the bands shifts the appearance of the Mott 
insulator phase towards deeper optical lattice potentials and stronger on-site interaction strength, 
as atomic coherence can build up via photon assisted tunneling in both bands. Varying the intra 
and interband on-site interactions leads to several different atomic phase configurations. There are 
even parameter regions where a mean field approach predicts concurrence of a Mott insulator state 
in one band, while atoms in the second band stay superfluid. 

PACS numbers: 37.10.Jk, 42.50.Pq, 37.30,+i 



I. INTRODUCTION 

Cavity Quantum Electrodynamics (CQED), i.e. the 
dynamics of an quantized electromagnetic field in a res- 
onator resonantly coupled to a system with quantized 
energy levels, is a central theoretical thought model to ex- 
hibit the physics of coupled quantum dynamics and mea- 
surements [1]. Experimental realizations of CQED were 
e.g. implemented for single Rydberg atoms in a super- 
conducting microwave resonator or for ground state al- 
kali atoms in high-Q optical resonators. In both cases the 
strong coupling limit could be reached, where light mat- 
ter coupling dominates the environmental decoherence. 
Later other implementations using solid state materials 
or ions followed ||. Recently, a Bose-Einstein Conden- 
sate (BEC) of ultracold atoms could be trapped within 
an optical high-Q cavity in the strong resonant coupling 
regime 0, 0, H[. As an important example, the role of 
cavity induced light forces on the atoms for the case of 
off-resonant coupling of a BEC in an optical cavity was 
also widely studied, e.g. in @, 0]. It seems just a matter 
of time that ultracold atoms in an optical lattice within 
a cavity will be achieved experimentally entering the full 
quantum many body domain of CQED 8] . 

In the present work we broaden our recent investiga- 
tions on excitons and polaritons in ultracold atom optical 
lattices, where we assumed a frozen Mott insulator state 
of the atoms in one or both bands Hence we study 

the motional quantum dynamics of a two-mode Bose gas 
in a prescribed optical lattice with resonant coupling of 
atomic excitations to cavity photons, to see under which 
conditions such a common Mott state will exist and be 
stable. In contrary to some recent studies on cavity medi- 
tated dipole forces Q , we assume resonance between the 
cavity mode and the atomic transition, so that cavity 
induced light forces are very small compared to the pre- 
scribed lattice potentials [llj and will thus be ignored 
in the following. In such cases collective electronic ex- 
citations (excitons) and cavity polaritons [H, [T3, EH EH 
depend on the atomic position distribution, which we 



have demonstrated studying defects in optical lattices 
[l4j . Here we will study the stability conditions of such 
a lattice itself, which was the basis of our previous cal- 
culations. While the optical lattice is formed by external 
classical fields far off resonance to any atomic transitions, 
the cavity field is represented by quantized single mode 
close to resonance to an atomic transition. In practise 
analogous coupling can be achieved by a resonant Ra- 
man transition between the two states, where one arm is 
connected to the cavity field [To] ]. 

In the following we concentrate on the modified phase 
diagram for such a system with particular emphasis on 
the influence of the cavity coupling on the transition from 
a superfluid to a Mott insulator phase. As we consider 
two boson modes for the ground and excited state atoms, 
quantum phase transitions can occur in each band sep- 
arately or in both bands together. Besides site to site 
hopping in each band an atom can also jump by changing 
bands hopping to another site and changing back to the 
original band. Hence even if the conditions for a Mott in- 
sulator are fulfilled in one band, the coupling can induce 
long range coherence in this band via indirect tunneling. 
Similar strong on-site repulsion between atoms in differ- 
ent bands could lead to an alternating state of atoms in 
being at different bands in neighboring sites. While for 
any concrete setup the hopping amplitudes and on-site 
interaction strengths can in principle be explicitly calcu- 
lated, that can be changed by modifying lattice depths 
and scattering lengths of the particles. Hence here we 
will keep them as free parameters to explore a wide pa- 
rameter space of such systems. While as first guess one 
might think of an optical cavity transition frequency for 
such a system, also implementation with microwave res- 
onators come to mind, where decay and decoherence are 
much smaller in general [lfjj . 

The paper is organized as follows. In section 2 we 
review the quantum phase transition for the single- 
component and two-components Bose-Hubbard model 
with interaction. In section 3 we add the cavity interac- 
tion between the two-components Bose-Hubbard model. 
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As above the corresponding phase-diagram is calculated 
in the mean field theory in section 4. A summary of the 
results appears in section 5. The appendix includes the 
derivation of the second order perturbation theory. 



II. ULTRACOLD ATOMS IN AN OPTICAL 
LATTICE 

A. Single-Component Bose-Hubbard Model 

Bose-condensation of atomic vapors and successive 
loading of the particles into an optical trap were ground 
breaking recent experimental achievements, demonstrat- 
ing unprecedented control of internal and external quan- 
tum state of particles. Cold atoms have been trapped in 
laser fields forming a standin g w ave which creates a peri- 
odic lattice of microtraps (l7l Il8l| . When loaded into such 
an optical lattice cold atoms still hop from site to site via 
tunneling. In addition, atoms occupying the same lattice 
site will repel each other due to collisional interactions. 
The stronger the laser field, the deeper the lattice and 
the slower the hopping rate of the atoms. At the same 
time the interaction energy between the atoms becomes 
stronger, as the atoms on one lattice site get more com- 
pressed. Thus laser light allows one to control the kinetic 
versus the repulsion energy of the quantum gas. In this 
way one can switch, e.g., a Bose gas from a weakly inter- 
acting gas, where atoms form a superfluid, to a strongly 
correlated regime, where atoms avoid each other and lo- 
calize into individual sites forming a so called Mott insu- 
lator phase. This cross over between two quantum phases 
at T = as a function of an external parameter is called 
a quantum phase transition, a phenomenon which in the 
meantime has been observed in many laboratories. 

Naturally this quantum phase transition from the su- 
perfluid to the Mott insulator phase was widely stud- 
ied both theoretically [l7], EH E(| and experimentally 
[3 120 . [2 ll| and it was found that the Bose-Hubbard 
model 22[ predicts most properties of this phase tran- 
sition very well. Trapping potentials can be experimen- 
tally controlled and using different initial BEC-densities 
the optical lattice can be filled with one, two or even 
more atoms per site on average by changing the relevant 
parameters. 

Mathematically, the single component Bose-Hubbard 
model is represented by 



H 



^E^ + Efe--") blh + ^bMbA. (1) 



The i summation is over the lattice sites, where b\ and 
bi are the creation and annihilation operators of an atom 
at site i, and £j is the atoms on-site energy, which in- 
cludes the external trap potential, /i is the chemical 
potential accounting for a possible exchange of atoms 
with an external atom reservoir. We assume here finite 
hopping (atom transfer) only between nearest neighbor 



sites, as indicated by the bracket, with the hopping pa- 
rameter J. Here U is the effective on-site atom-atom 
interaction strength, which is taken to be a repulsive and 
is characterized by the s-wave scattering length. Close 
to T = the atoms stay localized in the lowest Bloch 
band of the optical lattice potential. Experimentally the 
relative magnitude of parameters J and U can be con- 
trolled by changing the laser intensity (potential depth). 
In the limit of J ^> U the atom hopping among the lat- 
tice sites dominates the dynamics, and the atomic states 
are spread out over the whole lattice, so that we get a 
superfluid phase. In the opposite limit J <C U , the on- 
site atom-atom interaction dominates and the atoms are 
localized on individual sites, so that we have a Mott in- 
sulator phase. Between the two limits we get a quantum 
phase transition, which is nicely ^predicted by the Bose- 
Hubbard model presented above [13] ■ The Mott insulator 
phase is characterized by a fixed integer number of atoms 
per site. 

In figure (1) we plot the corresponding phase diagram 
for the plane (p/zJ) — (U/zJ) which is scaled by zJ, 
where z is the number of nearest neighbors. The super- 
fluid (SF) phase appears outside the three Mott-insulator 
(MI) phase regions, which are for one, two, and three 
atoms per site, that is n = 1,2,3. Here we assumed 
£j = 0. Such a phase diagram can be conveniently calcu- 
lated in a mean- field approach [HI, [24[ . 
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FIG. 1: Phase diagram: (fi/zj) vs. (U/zJ). The superfluid 
SF and the three Mott insulator MI regions, for n = 1,2,3, 
are shown. 



B. Two-Components Bose-Hubbard Model 

Let us extend this model to include two different inter- 
nal states of the atom, called ground and excited state. 
Hence we consider a boson gas of two-level atoms de- 
scribed by two kinds of bosons differing by their internal 
state. Equivalently also to different bands for atoms in 
the same internal state could be used. As before, the 
atoms are cooled and loaded on an optical lattice. In 



general the light forces of the lattice lasers are acting dif- 
ferently on ground and excited state atoms, and we get 
different optical lattice potentials for them as shown in 
figure (2). Quite generally we can, nevertheless, assume 
optical lattice potentials with minimums at the same po- 
sitions. For certain wavelengths, called magic laser wave 
lengths [25| , even identical potentials can be found, which 
has particular interest for building lattice atom clocks. 
For cold enough particles the ground and excited state 
atoms can be assumed to be localized in their first Bloch 
band. In this case the system is described by a two- 
component Bose- Hubbard model [26|, [2?j], given by the 
Hamiltonian 

H = -J a J2 b i b j- J °I2 c l c i 

(»,i> (»ij> 

+ t E b t b M + y E 44** + u eg e b\ c \b lCl 

i i i 

+ - ^9) b % + " A*e) 4«> ( 2 ) 

i i 

where b\ and 6,; arc the creation and annihilation opera- 
tors of a ground state atom at site i, respectively, with 
the on-site energy ef , and chemical potential fi g ; cj and Cj 
are the creation and annihilation operators of an excited 
state atom at site i, respectively, with on-site energy e|, 
and chemical potential fi e ; where e\ — ef + huj a , with 
uj a being the effective atomic transition frequency, which 
can include AC Stark-shifts relative to the free atomic 
transition frequency Here J g and J e are the hopping 
parameters of ground and excited state atoms, respec- 
tively. U g and U e are the on-site ground and excited 
state atom-atom interactions, respectively, and U eg is the 
on-site ground-excited atom interaction related to scat- 
tering between ground and excited state atoms. This 
extended model leads to a much richer phase diagram 
including also boundaries for the superfluid to the Mott 
insulator quantum phase transition. The Mott insulator 
phase with a fixed atom number per site can be of differ- 
ent form from above, i.e. we have one atom per site but 
it does not matter, if it is ground or excited state atom. 
Also it is possible to get phases with two atoms per site, 
with a number of options: two ground state atoms, two 
excited state atoms, or one ground and the other excited 
state atoms per site, and so on to higher atom number 
per site. 

Let us show some typical examples now. In figures (3- 
4) we plot the phase diagram for the plane (/i g /zJ g ) vs. 
(Ug/zJg), and (fx e /zJ e ) vs. {U e /zJ e ), which are scaled by 
zJ g and zJ e . In figure (3) we used for the scaled atom- 
atom coupling the number U eg /zJ g = U eg /zJ e = 15. 
The plot show the superfluid phase (SF) and the Mott 
insulator phase (MI) for one ground or excited state atom 
per site, that is n g = n e = 1, where the transition lines 
for ground and excited atoms are coincide. Figure (4) is 
with different atom-atom coupling, where U eg /zJ g = 15 
and U eg /zJ e — 20. Here the transition line (Le) is for 
excited state atoms, and (Lg) for ground state atoms. 
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Ground state optical lattice 



FIG. 2: Schematic plot of the ground and excited state optical 
lattice potentials. 

Beside the superfluid phase (SF), we get three other re- 
gions. The Mott insulator (MI) region with one ground 
state and one excited state atom per site. The (SM) 
region, where the excited state atoms are in the Mott 
insulator phase with one atom per site, and the ground 
state atoms are in the superfluid phase. The (MS) re- 
gion, where the ground state atoms are in the Mott in- 
sulator phase with one atom per site, and the excited 
state atoms are in the superfluid phase. In the previ- 
ous two cases of figures (3-4) we assumed ef = e\ = 0. 
In figure (5) we used for the scaled atom-atom coupling 
the number U eg /zJ g — U eg /zJ e — 15, with ef/zJ g = 
and ef/zJ e = 100. Here the excited and ground phase 
transition lines of figure (3) split. The phase diagram is 
calculated here by using the mean-field approach [26| . 




FIG. 3: Phase diagram: (n g /zJ g )vs. (U g /zJ g ), and (fj, e /zj e ) 
vs. (U e /zJ e ). The superfluid SF and the Mott insulator MI 
regions, for n g — n e — 1, are shown. We used U eg /zJ g = 

Ueg/zJe. = 15. 



III. AN OPTICAL LATTICE WITHIN A 
CAVITY 

In a next step we will now add resonant interaction 
with a cavity mode. Note that trapped ultracold atoms 
within a cavity were already experimentally studied by 
a number of experimental groups 0, 0, S EB| ; but as far 
as we know no optical lattices within a cavity were real- 
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FIG. 4: Phase diagram: (fi g /zJ g ) vs. (U g /zJ g ), and (fi e /zj e ) 
vs. (Ue/zJe). The SF and MI regions, for n g — n e — 1, are 
shown. The SM region is for the Mott insulator of excited 
atoms and superfiuid of ground atoms, and the MS region 
is for the Mott insulator of ground atoms and superfiuid of 
excited atoms. The full-line Lg is for the transition line of 
ground atoms, and the dashed-line Le is for the transition line 
of excited atoms. We used U eg /zJ g = 15 and U eg /zJ e = 20. 



rotating wave approximation, is 

H = J2(fi4bia + f* b\cHj, (3) 

i 

where fi is the coupling parameter, which is taken here 
to be of the electric dipole interaction. The first term 
represents the excitation of an atom from the ground 
state into the excited state by the absorption of a photon, 
and the second is for the jump of an excited state atom 
into the ground state by the emission of a photon. 




FIG. 6: Ultracold atoms in an optical lattice within a cavity. 
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FIG. 5: Phase diagram: {p, g j zj g ) vs. (U g /zJ g ), and (fi e /zj e ) 
vs. (Ue/zJe). The SF, MI, SM, and MS regions, for n g = 

15, with 



n e = 1, are shown. We used U eg /zJ g 
£ 9 JzJ g = and e e JzJ e = 100. 



U e g / zj e 



ized yet. Here we investigate a situation as depicted in 
figure (6), where we start from the previous system of op- 
tical lattice with two state atoms, which we now placed 
within a cavity [||. We assume that the atomic transi- 
tion is close to resonance with only a single cavity mode 
equivalently coupled to all atoms. Hence either mode 
function is constant over the lattice or we consider very 
long wavelength transitions (microwaves) [l(| , where the 
atoms are within a single wavelength. The mode Hamil- 
tonian is H c = e c a^a, where a^ and a are the creation 
and annihilation operators of a cavity mode of energy e c , 



respectively, and where e c 



Eg. The coupling be- 



tween the cavity mode and the atomic transition, in the 



As central goal of the present work we check how the 
phase diagram is changed as a result of the cavity induced 
interactions. Here important differences are expected due 
to the fact that the absorption and emission of cavity 
photons couple the two kinds of bosons and allow for 
nonlocal energy transfer between the atoms. A transition 
from an excited to a lower state atom at one site with 
photon creation can be compensated by an excitation 
at any other site by photon absorption. In principle this 
process can be mediated by a virtual photon, which never 
actually appears in the mode. 

In particular we examine, how the superfiuid to the 
Mott insulator phase transition is affected by such cou- 
pling to the cavity photons, which should allow to build 
long range correlations very effectively. At this point our 
approach here is mainly analytical based on mean field 
theory for atoms and the field. We apply method as it 
appears in [23j, which was successfully used in [24j to 
one component Bose-Hubbard model, and in [26[ for two 
components case, where the superfiuid and the Mott in- 
sulator phases predicted to a good approximation, and 
with their quantum phase transition. We assume long 
life times for both the atomic excitations and the cavity 
photons, where their line widths are taken to be smaller 
than the excitation-photon coupling. Hence in the follow- 
ing calculations we neglect the damping for the excita- 
tions and the photons. Furthermore, we consider a cavity 
with a fixed number of photons, which can be achieved 
by the combination of an external pump and the photon 
damping. 

A two-components Bose-Hubbard model within a cav- 
ity in the resonant regime is represented by the Hamilto- 
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and 



(id) (id) 
i i 
i i 

+ u ea « + E (/ ^ a + f* at b \ c ) . ( 4 ) 



where here \x g — p, g — e 9 and /i e = p, e — e e , with /2 ff 
and /i e are the ground and excited state atoms chemical 
potentials, respectively, ef and e\ are taken to be site 
independent, hence we dropped their index i. At pump- 
cavity mode resonance the chemical potential for cavity 
photons is zero, that is p, c = 0. Furthermore, the photon- 
excitation coupling / is assumed to be site independent, 
which is the case for homogeneous optical lattice parallel 
to the cavity mirrors. We also defined the mean atoms 
and photons number operator, by nf = b\b{, n\ = cjc,-, 
and n c = a^a. 



IV. PHASE DIAGRAMS IN THE MEAN FIELD 
THEORY 

In the mean-field theory presented in [23| the atoms 
are subject to the mean-field of the neighboring sites and 
of the other kind of bosons. Hence the Hamiltonian can 
be separated into on-site terms and the rest, which can 
be written as a sum of a free part and an interacting part. 
Then the interaction part can be treated in perturbation 
theory. 

In the mean-field theory we have 

(h) = (b\) = J^=4, g) (a) = (4) = V< = fa. (5) 



We use the consistent mean-field to obtain the hopping 
and the excitation-photon coupling terms, by 



btbi = (b\ + bi) 4> g -^ g , 



c]a = (ct + c, 

c]h = hi (b e + c: 



(6) 



The Hamiltonian is now given by H = Hq + Hj , where 
H o = -Mg E n i ~ ^ e E n i + £ c nC + U eg E 



n\n\ 



^E n f( B ?- 1 )+TE n ' , w- 1 )' (7) 



Hi = -zJ g J2{ b t + bi ) ^ + zJ 9j2' 

i i 

- zj e ( c| + cl) 4>e + z j e 4>l 

i 

+ E f a ( 6 « 

i 



The Hamiltonian is on-site separated, where Hq 
E^o, with 

Hq = -fj, g n g - lJL e n e + £ c n c + U eg n g n e 



U, 



U, 



+ -f ng K - 1) + ~f ri e (n e - 1) . (9) 

We formally defined n c = n c /N in dividing by the num- 
ber of lattice sites N. Hence we dropped the site index, 
and to write n g = nf, n e = n|, and n c = ff . 

The interacting part is Hj = J2iH}, with H\ = H l hop + 
H* av , where, after dropping the index i, we get 

HL P = ~zJg {b ] +b) 4> g + zJ g 4> 2 g 

- zJ e (c f +c) (/>e + zJ e (10) 



and 



HI 



= f a (b <f> e + c+ (j> g - (f> g (t> e ) 

+ f* 0e + C0 fl -^e)- (11) 



Now we use perturbation theory, where in the mean- 
field theory the hopping and the excitation-photon cou- 
pling is a perturbation. The system ground state is repre- 
sented in occupation number states by |0) = \n g ,n e , n c ). 
To the zero order we get 

(0|i?o|0) = i?o 0) = -fig n g - iJt e n e + e c n c + U eg n g n e 



U, 



^ n g {n g - 1) + ^ n e (n e - 1) . (12) 

At ground state we have dE^ /dn e = BEq /dn g = 0, 
which gives 



He = U e [n 



l 



(13) 



fig = U g l^rig - - j + Ueg n e , 

which are solved to get 

U g (2 Me + U e ) - U eg (2/Xg + U g ) 

2 (UeUg - Utg) 

= Ue(2Vg + Ug)-Ueg (2fjL e + Ue) 

2{U e U g -U%) ■ 1 ' 
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Here n g and n e need to be integers, hence we define n g = 
n g + a and n e — nP e + (3, where —1/2 < (a, [3) < 1/2, or 
(n°-l) < (n s -l/2) < n° g and (n° e -l) < (n e -l/2) < n° e , 
which is solved for fi g and fi e , and which insure n g and n e 
to be integers. If n g or n e is zero, then the ground state 
is unstable due to the interaction between different atom 
kinds, U eg . At n g — or n Q e = 0, for the inequality to hold 
we need the conditions U g > 0, U e > 0, and U e U g > U 2 g 
[28| . Therefore, to avoid instability we assume ground 
state with non-zero and positive atom number, that is 
(n g ,n g ) > 0, or to limit the discussion to the previous 
stability conditions. 



In a second step we are now including higher order 
corrections to the ground state energy induced by the 
hopping among nearest neighbor sites and by the excita- 
tion coupling to the cavity photons, in calculating higher 
order terms of the perturbation series. As the first or- 
der term is zero Eq = 0, we thus have to calculate the 
second order term 



(2) 



K0|ifz|n)f 



(0) 



(15) 



The details of the calculation appear in Appendix A, 
which yield 



n g + 1 



fig - Ug Tig - U eg Tl e 



-fig + Ug (Tig + U e g 



l/l 2 



(n e + 1) n c 



fle + £c - U e n e - U eg 7l g 

n e (n c + 1) 

-[ie - £ C + U e (n e - 1) + U eg n g 



zJ e + Z J. 



1 T 2 



n e + 1 



fJ*e U e TL e U e g Tig 



+ 



in 



fle + Ue (n & ~ 1) + U eg Ug 
UnTle 



-fig + e c + Ug (n g - 1) + Ueg n e 



(n g + 1) (n c + 1) 

llg-£ C - Ug Tig - U eg Tl t 



€ (16) 



According to Landau theory, at phase transition the fac- 



tors of (j)g and 4>2 equal to zero, that is 



1 + z.J a 



n 9 + 1 



(~^9 ^9 ^9 ^ e 9 ^ e 

Hi 

-Ug + U g (rig - 1) + Ueg r 



+ 



\fl 

Zj a 



(lie + l)n c 



H e ~\~ S C U e 7l e U e g Tig 

n e (n c + 1) 

-fl e ~ £c + U e (n e - 1) + U eg Tig 

n„ + I 



0. 



1 + Zj e 



f-^e U G Tl e U e g Tig 

ne 

-fle + U e (n e - 1) + U eg n g 

l/l 2 \( n g n c 

ZJe 



-fig +£ C + U g (n g - 1) + U E g n e 



(n g + 1) (n c + 1) 

v fl>g £c Ug Tl g U eg Tl e 

In the limit zJ g , zJ e 3> |/| 2 we get 



1 + zJ„ 



(17) 



1 



Mfl ^9 ^9 ^ G 9 ^ e 

Kg 

-fig + Ug (Ug - 1) + U eg T 



1 + Z,J e 



n e + l 



fle -U e Tl e - U eg Tl, 
II, 



, ( 18 ) 
-fl e + U e (n e - 1) + Ueg Ug / 

which is the result for two-component Bose-Hubbard 
model in neglected the cavity effect, and which leads to 
the results in figures (3-5). 

In the limit where the excitation-photon coupling is 
much larger than the hopping, that is zJ g ,zJ e <C |/| 2 , 
we get 

|/| 2 IV K + iK 



+ 



i , { 

Zjg \^fi e ~\~ £ c ' U e Ti e ■ U e g Tig 

n e (n c + 1) 

-fle -£ C + U e (n e - 1) + Ueg U g 

1+ m 2 



0. 



-fig +e c + Ug (Tig ~ l) + U eg 



(n g + 1) (n c + 1) 



0. 



(19) 



v fLg £c Ug Tig Ueg Tie 

which we are going to examine rigorously. We solve for 
Jig = fig + e 9 and fl e = fi e + £ e , to get the rescaled results, 
which are rescaled in dividing by zJ g and zJ e , and after 
dropping the bars, we get 



fig 



fit 



^- = (sg+el)+L g ± 

ZJn 



- ^g, 



t^- = (i e -sl) + L e ± 

ZJ e £ 



AK e ,(20) 
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where 



L„ = 



with 



and 



1 



F 



jf a = U g [n g - - ) + U% n e - - (n g + n c + 1) , 
-\)+Ut g n 9 --(j, - ,,, : (21) 



G„ 



G e = 



if, = 



G^ 

Zjg 

Ge_ 

zJ K 



U g + F (n g + n c + 1) , 

f) e + F (n c - n e ) , (22) 

= FU- g (n g + 1) (n c + 1) , 



K K 



(zJ e 



= FU e (n e + l)n c , 



(23) 



phase to large atom-atom interactions, which is here one 
order of magnitude larger. Also we get that the ground 
state transition line is shifted relative to the excited state 
one. 

In figures (8-11) we plot the scaled chemical potentials 
/ig and [i e as a function of different parameters, for n g = 
n e = 1. In figure (8) the plot is as a function of U eg , 
in using e c = 100, F = 25, U g = U e = 250, and n c — 
1. In figure (9) the plot is as a function of e c , in using 
U 



eg 



15, U g = U e = 250, F = 25, and n c = 1. In figure 

(10) the plot is as a function of F, in using U eg = 15, 
U g = U e = 250, e c = 100, and n c = 1. It is clear 
from figure (10) that the Mott insulator phase appears 
only for a limited range of coupling parameters. In figure 

(11) the plot is as a function of n c , in using U eg = 15, 
U g = U e = 250, e c = 100, and F = 25. Also here the 
Mott insulator phase is obtained for a limited range of 
mean cavity photon number per site, which ranges from 
few cavity photons up to about one photon per lattice 
site. 



where 



l/l 5 



(24) 



with U g = %,U e = jfe, and U ig = Ut g = % 



also ga c = 7£> 



= , and we have e„ = — f - , e e 

ZJ C ' 9 zj„ ' c 



As we concentrate in the following on the influence of 
the coupling to cavity photons, we assume J g = J e = J, 
hence we drop the tildes and all the parameters to be 
considered as scaled in dividing by z J. Furthermore, we 
assume zero ground state energy, e g — 0, and the cavity 
mode is in resonance to the electronic excitation, that is 

In figure (7) we plot the phase diagram for the plane fi g 
vs. U g , and the plane \x e vs. U e . We used for the scaled 
atom-atom coupling the number U eg — 15. The plot 
shows the superfluid phase (SF) and the Mott insulator 
phase (MI) for one ground or excited state atom per site, 
that is n g = n e = 1, and where we used for the cavity 
photons n c = 1. The transition line (Le) is for excited 
state atoms, and (Lg) for ground state atoms. Beside 
the superfluid phase (SF) and the Mott insulator (MI) 
regions, we have the (SM) region, where the excited state 
atoms are in the Mott insulator phase with one atom per 
site, and the ground state atoms arc in the superfluid 
phase. The (MS) region, where the ground state atoms 
are in the Mott insulator phase with one atom per site, 
and the excited state atoms are in the superfluid phase. 
In the figure we assumed e c — 100 and for the excitation- 
photon scaled coupling we used F — 25. We compare the 
present results of figure (7) to that of figure (5) for the 
phase diagram of two-components Bose-Hubbard model 
without a coupling to cavity photons. We deduce that 
the effect of the coupling to cavity photons, with one 
photon per lattice site n c = 1, and for strong excitation- 
photon coupling of F = 25, is to shift the Mott insulator 
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FIG. 7: Scaled phase diagram: 
for n g = n e = 1. We have e c = 100, F = 25, U eg = 15, 
n c = 1. The dashed line is for the excited state atoms, 
the full line is for the ground state atoms. 



Ue, 

and 
and 



In figure (12), beside the lines for the ground and ex- 
cited state transitions from the superfluid into the Mott 
insulator with one atom per site (the the upper full line 
is for the ground state atoms, and the lower full line 
is for the excited state atoms), we plot the transition 
line for the case of two ground state atoms per site (the 
upper dashed line), when the excited state atoms are 
with zero atoms per site. The plot is for fi g vs. U g , 
and /it, e vs. U e . We used the parameters: e c = 200, 
F = 25, and U eg = 50, where the plot is for the two 
cases of n g = n e = 1 and n g = 2, n e = 0. Here 
we assumed n c = 1. At these parameters the Mott 
insulator for n g = n e = 1 is separated from that for 
n g = 2, n e = 0, as no crossing appears between the full 
and the dashed lines. In figure (13) the plot is for the 
same parameters but now with U eg = 500. It is clear 
that for large U eg the transition lines of n g = n, 



and n a 



n e = 1 

are crossed, and we get correlations 
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FIG. 8: Scaled phase diagram: fi g and /i e vs. U eg , for n g — 
n e = 1. We have e c = 100, F = 25, t/ 9 = t/ e = 250, and 
n c = 1. The dashed lines are for the excited state atoms, and 
the full lines are for the ground state atoms. 
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FIG. 10: Scaled phase diagram: (/i g ,/i e ) vs. F, for n g — n e = 
1. We have U eg = 15, U g = U e = 250, e c = 100, and n c = 1. 
The dashed line is for the excited state atoms, and the full 
line is for the ground state atoms. 
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n c = 1. The dashed lines are for the excited state atoms, and 
the full lines are for the ground state atoms. 



between the Mott insulator with one ground and one ex- 
cited atom per site and that with two ground and zero 
excited atoms per site. Figures (14-15) are for the tran- 
sitions with n g — n e — 1 and n g = 1, n e — 0. In figure 
(14) we used U eg = 500, and in figure (15) U eg = 50, 
where the other parameters are the same as before. For 
small U eg we get strong correlations between the Mott in- 
sulator phase with n g = n e = 1 and n g = 1, n e = 0. For 
larger U eg the correlation appears only for large U g and 
U e , but for small ones the transition lines are separated. 
In figure (16) the plot is with U eg — 250. Here the two 
full lines are for the transitions with n g = n e = 1, the 
lower dashed line is for n g = 1, n e = 0, and the upper 
dashed line is for n g = 2, n e = 0. We conclude that no 
crossing between the transition lines of n g = 1, n e = 
and n g = 2, n e = is obtained. Furthermore, in fig- 
ure (17) we plot the phase diagram for the three tran- 



sitions with (n g 



1, n e 



0), and 



FIG. 11: Scaled phase diagram: (/i g ,/i e ) vs. 



1. We have U e 



15, U g = U e = 250, e c = 100, and 



F = 25. The dashed line is for the excited state atoms, and 
the full line is for the ground state atoms. 



the transition with n g = n e = 1, where the upper for the 
ground state atoms and the lower for the excited state 
atoms. The upper dashed line for the ground state atoms 
with n„ = 1, n e = 0, and the lower dashed line for the 



= 1. Here we used 
Each zone in the 



excited state atoms with n g = 0, n e - 
also e c = 100, F = 25, and U eg = 30. 
diagram can be explained as before, just to note here the 
correlations between the different Mott insulator regions. 
For illustration in figure (18), for the same parameters, 
we plot the phase diagram for the three transitions with 
(n g = n e = 1), (n g = 2, n e = 0), and (n g = 0, n e = 2). 
The two full lines are for the transition with n g = n e = 1 . 
The upper dashed line for the ground state atoms with 
n g = 2, n e = 0, and the lower dashed line for the excited 



(n g = 0, n e = 1). As before the two full lines are for state atoms with n g = 0, 
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FIG. 12: Scaled phase diagram: /i 9 vs. U g , and /i e vs. U e , for 
n g = n e = 1 and n g — 2, n e — 0. We have e c = 200, F — 25, 
(7 e9 = 50, and n c = 1. The dashed line is for the transition 
with n 9 = 2, n e = 0, and the full lines are for the transition 
with n„ — n e = 1. 



FIG. 14: Scaled phase diagram: /i 9 vs. [7 9 , and /i e vs. {/ e , for 
w 9 = n e = 1 and n g = 1, n e = 0. We have e c = 200, F — 25, 
(7 e9 = 500, and n c = 1. The dashed line is for the transition 



with n 9 = 1, 
with n 9 = n e 



n e = 0, and the full lines are for the transition 
= 1. 
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FIG. 13: Scaled phase diagram: fi g vs. U g , and fi e vs. U e , for 
rig = n e = 1 and n 9 = 2, n e = 0. We have e c = 200, F = 25, 
I7 e9 = 500, and n c = 1. The dashed line is for the transition 
with n 9 = 2, n e = 0, and the full lines are for the transition with n 9 

with n„ 



FIG. 15: Scaled phase diagram 

n, 



fj, g vs. U g , and /i e vs. Ue, for 
n e = 1 and n g = 1, n e = 0. We have £ c = 200, F = 25, 
(7 eo = 50, and n c = 1. The dashed line is for the transition 
1, n e = 0, and the full lines are for the transition 
n e = l. 



V. SUMMARY 

We calculated the quantum phase diagram of a dilute 
gas of ultracold atoms occupying two bands of an opti- 
cal lattice within a cavity. The atoms in the ground and 
excited internal states are treated as two kinds of bosons 
modeled by a two-component Bose-Hubbard Hamilto- 
nian. Using a well established mean field based approach 
we analytically calculate and plot the parameter condi- 
tions for the transition from the superfluid to the Mott 
insulator phase for a wide range of parameters. When 
the optical lattice is located within a cavity with a single 
cavity mode close to resonance with the transition be- 
tween the two internal electronic states, the possibilities 
for tunneling of the atoms increase. Hence, recalculating 
the phase diagram in the strong coupling regime we found 
that the effect of coupling to cavity photons indeed shifts 



the Mott insulator region of the phase diagram towards 
much larger atom-atom interaction parameters. This re- 
sult is due to the fact that cavity photons tend to delo- 
calize the excited atoms in the optical lattice. Besides 
the phase transition borderline for one atom per site, a 
more general approach allows to obtain transition lines 
for higher numbers of ground and excited state atoms 
per site showing a rather complex and overlapping gen- 
eral phase diagram. 

Nevertheless, for strong enough interactions despite all 
these couplings, a parameter region of a Mott insulator in 
both bands, where the atoms are not moving can be iden- 
tified. Such a state was assumed as basis in our previous 
calculations on excitons and cavity polaritons in optical 
lattices. These models now could also be extended to 
mixed lattice phases. 

The results of the present paper should guide future ex- 
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FIG. 16: Scaled phase diagram: /i 9 vs. U g , and /j, e vs. U e , for 
(n 9 = n e = 1), (n g = 1, n e = 0), and (n 9 = 2, n e = 0). We 
have e c = 200, F = 25, (7 e9 = 250, and n c = 1. The upper 
dashed line is for the transition with n g — 2, n e = 0, the 
lower dashed line is for the transition with n g — 1, n e = 0, 
and the full lines are for the transition with n g — n e = 1. 
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FIG. 17: Scaled phase diagram: /i 9 vs. i7 9 , and vs. U e , 
for (n 9 = n e = 1), (n 9 = 1, n e — 0), and (n 9 = 0, n e = 1). 
We have e c = 100, F = 25, U eg = 30, and n c = 1. The upper 
dashed line is for the transition with n g = 1, n e = 0, the 
lower dashed line is for the transition with n 9 = 0, n e = 1, 
and the full lines are for the transition with n g = n e = 1. 



periments on an optical lattice within a cavity as they are 
currently set up in a number of laboratories worldwide. 
We give a qualitative overview and even first quantita- 
tive estimates for required system parameters in order 
to achieve a quantum phase transition towards a per- 
fectly ordered state. Quantitatively our results can be 
improved in calculating higher order terms of the pertur- 
bation theory and in going beyond the mean field theory. 
Interestingly some preliminary numerical calculations on 
small system sizes show already completely different pos- 
sibilities of atomic order like e.g. alternating occupation 
probabilities between neighboring sites. These phases 
cannot be captured by our current approach and thus 
leave much room for future work. Similarly the inclusion 
of atomic dipole-dipole interactions in form of excitons 



1000 



800 




PoO 200 ^OOj 400 500 

9' e 



FIG. 18: Scaled phase diagram: /i 9 vs. U g , and fi e vs. U e , 
for (n 9 = n e = 1), (n g = 2, n e = 0), and (n 9 =0, n e = 2). 
We have e c = 100, F = 25, U eg = 30, and n c = 1. The upper 
dashed line is for the transition with n g — 2, n e = 0, the 
lower dashed line is for the transition with n 9 =0, n e = 2, 
and the full lines are for the transition with n g = n e — 1. 



will also lead to new physics beyond the Bose Hubbard 
model. 
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APPENDIX A: THE SECOND ORDER 
PERTURBATION THEORY 

Here we calculate the second order term of the 
perturbation series Eq. p^|) . We start to calcu- 
late the matrix elements (0\Hi\n), namely we want 
to calculate the matrix elements of the operators: 
b, b\ c, c\ a, cjT, ab, a^c, ac* , ofib\ where |0) = 
\n g ,n e ,n c ). They have non-vanishing matrix elements 
between the states 



11) = 


\n g + l,n e , 


n c ) 


, 12) = 


n g - l,n e ,n c ), 


|3) = 


\n g ,n e + 1, 


n c ) 


, K) = 


n g ,n e - l,n c ), 


|5) = 


\n g ,n e ,n c - 


-1) 


. 16) = 


n g ,n ei n c - 1), 


|7> = K 


+ l,n e ,n c - 


-1) 


. |8> = 


n g ,n e + l,n c - 1) 


19} - K, 


n e - 1, n c - 


-1) 


■ [10) = 


\n g - l,n e ,n c - 1 



(Al) 
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with the energy differences 



and 



(0) 
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(o) 



(0) 


(0) 


(0) 


(0) 



25<°> = 
£f = 
M 0) = 



E 



(0) 
4 

(0) 



Ek> = 



E, 



(o) 
o 

(0) 


(0) 


(0) 



(0) _ 

6 _ 

4 0) = 

E^ = 



where -E'q '' is defined in Eq. (fT2|) . The matrix elements 
are 



(0) 
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(0) 
10 



Ms 9 ^9 eg ^g.-) 

-fig + Ug (Ug - 1) + U e g H e , 
jl e U e 7l e U e g Tig-, 

-fi e + U e (n e - 1) + U eg n g , 
-e c , 

jig Eq Ug Tig Ueg TLe-> 
fj, e ~\~ S c U e Tl e Ueg Tig-, 

-fi e - e c + U e (n e - 1) + U eg n g , 

-Mfl +£ C + Ug (fig - 1) + U eg U e , 



(A2) 



n g - 


h 1 


'n e - 


h 1 



(0|6t|2) = V%, 



(0|o|5) = y/n c + 1 
(0|6a|7) = + lVn-c + 1 



| C t|4) = 
| t| 6 ) = 

I™tl8\ - 



91 



|8) = V^T+T 



lac 1 



rV^TTT , (0|aV|10) = 



(A3) 



= Iff 



(n e + 1) n c 



fle + £c - U e n e - U eg 7l g 

n g n c 

-Ms + £ c + U g (n g - 1) + U eg n e 
(ng + 1) (n c + 1) \ ^ 2 

fig - E C — Ug Ug ~ U e g Tl e J " 

n e (n c + 1) 

-\i e - e c + U e (n e - 1) + U eg Ug 

c i2 ±2 ( n c + 1) ,2/2 
-VeVg " §Ag 



(A5) 



The second order correction for the energy are 

n„ + 1 



E {2) - z 2 7 2 6 2 

a hop - z J g Vg 



fig - Ug 7lg ~ U e g 71, 



+ z z J, 



fig + Ug (n g - 1) + Ueg n t 

He + 1 



2 r2 i2 



fl e - U e n e 



Ueg rig 



-fl e + U e (n e - 1) + U eg Tl c 



zJ. 



9 Yg 



zJ e 



(A4) 



The terms of the order (j> 2 (f> 2 can be neglected at the 
second order. 
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